The effect of limiting resources in aging populations 
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Abstract: The concept of a carrying capacity is essential in most models to prevent unlimited growth. Despite 
the large amount of deaths it introduces, the actual influence of the Verhulst term in simulations is often times not 
accounted for. Generally, it is treated merely as a scaling parameter that functions to keep simulated populations 
within computer limits. Here, we compare two different implementations of the concept in the Penna model - 
Vehulst applied to all individuals (VA) and to newborns only (VB). We observe variations in certain model features 
when random Verhulst deaths are restricted to a single age group. 
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1 Introduction 

The Penna model [1] is a popular tool for simulat- 
ing aging in biological populations [2]. Since its in- 
troduction in 1995, more than 200 papers have been 
written utilizing the model Q. Much of its suc- 
cess is attributed to its simplicity and ability to re- 
produce universal features of much more complicated, 
real phenomena [4] such as Gompertzian mortality 
JU [5), catastrophic senescence of the Pacific salmon 
(6), sexual reproduction 0, speciation HOES, knowl- 
edge [10], and so forth. The model is based on the mu- 
tation accumulation theory of senescence which states 
that the strength of natural selection declines with age 
and a population can accumulate harmful mutations 
that only have late life effects [11]. That is, individu- 
als die from bad genes which reveal or express them- 
selves at later stages in life. 

The effects of finite resources are usually taken 
into account by setting a maximum sustainable popu- 
lation size - the environmental carrying capacity, K. 
It is generally believed that organisms compete for 
available food, space, and other necessities in order 
to stay alive. Thus, in computer simulations, the car- 
rying capacity concept usually takes the form of a sur- 
vival probability, also known as the Verhulst factor, 
V = 1- N(t) /K. Here, N(t) is the size of the pop- 
ulation at time t. The Verhulst factor has been used to 
demonstrate different phenomena such as the Eve ef- 
fect [12] and coevolution [131. But its main function 



is to limit growth lPT4l . 

In the original Penna model HMOl, the Verhulst 
term kills at random, regardless of fitness and age - 
VA implementation. Because random deaths in nature 
hardly play a significant role in population dynamics, 
this approach is not very well justified lfl6l . An alter- 
native was later introduced applying the Verhulst fac- 
tor to newborns only (VB implementation) lfT7l . This 
was done to avoid the accidental killing of healthy 
individuals lfl8l . Furthermore, since majority of the 
population consists of newborns, limiting its size is 
enough to prevent exponential growth. 

Whether VA or VB, a greater part of the popula- 
tion is killed because of the Verhulst factor. The ef- 
fectiveness of invoking a finite carrying capacity in 
capping simulated populations is undeniable. How- 
ever, extra caution must be made so as not to confuse 
Verhulst effects with particular model features; other- 
wise, results can be misleading [ 18 ]. 

This work investigates the consequences of vary- 
ing Verhulst influence on the average characteristics 
of simulated populations. In particular, we examine 
the age demographics and survival rates of the steady 
state populations resulting from VA and VB imple- 
mentations of the Penna model. 



2 The model 

The Penna model incorporates age structure in con- 
ventional population modeling via bit-handling tech- 
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Figure 1: (A) Mean population at steady state and (B) 
the corresponding standard deviations relative to the 
initial size (No = 20, 000) of populations obtained 
from a VA implementation with T = 3 and R = 8. 



niques [1]. Here, individual characteristics are stored 
in a string of binary numbers (genome) that is defined 
upon birth. The string contains genetic information 
- zeroes are healthy genes and ones are bad genes or 
diseases. It is then read in sequence at a rate of one 
bit per time interval. This is done at the beginning of 
each iteration. Whenever a new bit is read, the indi- 
vidual's age is increased by one. A one on the ith bit 
means that the individual will experience the effects 
of a disease starting at age i. 

An individual suffers a genetic death when the 
sum of active harmful mutations reaches the threshold 
value, T. The death age is determined for each indi- 
vidual at birth. When the total number of bad muta- 
tions in the genome is less than T, the individual dies 
at age equivalent to the bit-string length, L. Hence, 
the length of the string dictates the maximum lifes- 
pan. The concept of finite resources introduce ran- 
dom deaths into the simulation. Thus, for some, actual 
death may happen at an earlier age due to the Verhulst 
factor. 

Surviving individuals reproduce only once, at age 
R. The number of offspring generated by each parent 
is given by the birth rate, B. The newborns copy the 
genes of the parent and acquire one additional muta- 
tion that is set at a random location. In the study of the 
influence of heredity and mutations on the evolution 
of a population, only those with deleterious effects are 
often times considered. This is because harmful mu- 
tations are many times more frequent in nature than 
beneficial ones |fl9l. 
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Figure 2: Time evolution of the population associated 
with B = 4, T = 3 and R = 8. Values are plot- 
ted as factors of the initial size, Nq. Note that in our 
simulations, K = 10 x Nq. 



Simulation details: Genetic information is 
stored in 32-bit long genomes (L = 32). We begin 
our simulations with 20,000 perfect newborns (no bad 
genes). For VA and VB runs, the value of the car- 
rying capacity is set ten times the initial population 
(K = 200, 000). At each iteration, age is updated first 
and deaths come immediately after. This procedure 
sets a maximum reproduction age, R rna x = L — 1, and 
for the VB case, restricts the action of the Verhulst fac- 
tor to those at age 1 . To minimize fluctuations, demo- 
graphic data presented in the next section are averaged 
over the last 300 iteration. At these times, simulated 
populations are already at equilibrium. 



3 Results and Discussion 

In Penna simulations, saturation time is influenced by 
the reproduction age. Those with higher R values take 
longer to equilibriate [ 1 ]. The mutation threshold and 
the birth rate, on the other hand, affect population size. 
Total population generally increases with organism's 
tolerance for deleterious mutations. Figure UK shows 
a rise in mean population with birth rate. For these 
VA runs, the values saturate just under twice the initial 
population size. The magnitude of fluctuations (given 
by the standard deviation in Fig.QJ}) increases linearly 
with birth rate. The plots also show a maximum B 
value that yields a nonzero steady state. Large fluctu- 
ations associated with high birth rates sometimes lead 
to extinction. A sudden rise in size, beyond the set car- 
rying capacity, results in the accidental killing of the 
entire population. The same is observed in VB sim- 
ulations. Fluctuations in the transient part cause pre- 
mature extinction even at relatively lower birth rates 
(Fig. [2]>. Note that although the effect of the Verhulst 




Figure 3: Behavior of the survival rates when B = 3, 
R = 8: (A) T = 3 and (B) T = 27. 

factor is restricted to a single age group, the strength 
of this external death term is not diminished. In fact, 
the number of random deaths is doubled in the VB 
case [18]. Thus, the range of B values that result in 
nonzero steady states is smaller for the VB implemen- 
tation. It is for this reason that later comparisons of 
VA and VB simulation runs use different B values. 

Recall from Section 2 that we restrict reproduc- 
tion to a single age R. This one-time breeding strategy 
is characteristic of semelparous organisms. The aging 
process associated with semelparity is generally catas- 
trophic lUEZl. For the case of reproduction occuring 
only once at R = 8, we plot the survival rates as a 
function of age (Fig. [3]). In aging studies, the survival 
rateQis related to organism fitness. The abrupt decline 
in fitness happens at age R+ 1, immediately after re- 
production, when T < R; otherwise, when T > R, 
the drastic loss in functional abilities is observed at 
the age equivalent to the threshold [20]. Notice fur- 
ther that the survival rates associated with the VA case 
in Fig. [3}3 appear to have a period of R. This behav- 
ior is not a general feature of the model, but rather a 



'the ratio between the number of individuals with age a at 
time t and the number of individuals with age a — 1 at time t—1, 
JV„(t)/JV„_i(t-l) 



Figure 4: Distribution of the population to its com- 
ponent ages, normalized about the average number of 
births. Here, R = 16: (A) T = 1 and (B) T = 20. 
The B values used are 10 and 3 for the VA and VB 
simulations, respectively. 

consequence of our choice of initialization and Ver- 
hulst implementation [20]. The dips correspond to 
times when the Verhulst killing effect is strongest, felt 
by ages present at times when the population size is 
largest. 

The Penna model described in fl] [151 follows 
Gompertz law (2j[5l which predicts an exponential in- 
crease in mortality with age. Consequently, the distri- 
bution of the population to its component ages shows 
an exponential decrease. Note that these populations 
are VA and breed repeatedly, starting at age R until 
death. For our current implementation, we obtain the 
same exponential decrease in the age demographics 
when the Verhulst factor is applied to all ages. How- 
ever, when random deaths imposed by the Verhulst 
factor is restricted to newborns only, we lose the Gom- 
pertzian structure (Fig. [4]). 

Figure [5] presents the demographic data associ- 
ated with Verhulst-free populations. These are re- 
sults of simulation runs that do not account for en- 
vironmental restrictions. That is, there are no random 
killings by the Verhulst factor and all deaths are due 
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Figure 5: Demographic statistics obtained for the 
Verhulst-free case: B = 1, T = 20 and R = 16. 

to genetic reasons. Age demographics and survival 
rates obtained for the Verhulst-free case show no var- 
tiation in mortality rate and physiological condition 
with age, respectively. These are two of the criteria 
identified by Caleb Finch [21] for negligible senes- 
cence. These describe the very slow aging process ob- 
served in coldwater fish, bivalves, turtles, whales, and, 
fairly recently, in the naked mole-rats |[22l . Note that 
within the framework of the Penna model, Verhulst- 
free demonstrations are limited to threshold values 
that are greater than the reproduction age (T > R) 
with one semelparous birth per adult (B = 1) [23]. 
So far, these are the only parameter sets that lead to 
nonzero steady states even without random deaths by 
the Verhulst factor. 



4 Summary and conclusion 

Population dynamics models provide quantitative es- 
timates of the population size at some time t, given 
its initial size, growth rate, and the carrying capac- 
ity ll24l . The carrying capacity incorporates all possi- 
ble interactions among individuals in an environment 
with limited resources into a single parameter. Its 
value determines the size of the steady state popula- 
tion [25]. Without the concept of a carrying capacity, 
populations behave exponentially. Thus, most models 
deem this parameter necessary to avoid the exponen- 
tial growth. 

In computer models, the concept generally takes 
the form of a Verhulst factor - a random death proba- 
bility. Despite the fact that majority of deaths in simu- 
lations is attributed to this term, its effect on the prop- 
erties of simulated populations is rarely considered. 
The choice of implementation of this random death 
strategy unfortunately has an unexpected impact on 
the population's genetic profile lfT8l . In addition to 



this, we find that restricting the action of the Verhulst 
term to newborns only (VB implementations) alters 
some of the basic model features which was demon- 
strated imposing the death probability to all individu- 
als regardless of fitness and age (VA implementation). 
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